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FRACTIONAL DERIVATIVES AND FRACTIONAL POWERS OF 

OPERATORS 
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Abstract. We give stability and consistency results for higher order Griinwald- 
type formulae used in the approximation of solutions to fractional-in-space par- 
tial differential equations. We use a new Carlson-type inequality for periodic 
Fourier multipliers to gain regularity and stability results. We then generalise 
the theory to the case where the first derivative operator is replaced by the 
generator of a bounded group on an arbitrary Banach space. 



1. Introduction 

In a series of articles Meerschaert, Scheffler and Tadjeran [T|)J H3 [TH1 HH1 HO] ex- 
plored consistency and stability for numerical schemes for fractional-in-space partial 
differential equations using a shifted Griinwald formula to approximate the frac- 
tional derivative. In particular, in [20] they showed consistency if the order of the 
spatial derivative is less or equal to 2. They obtained a specific error term expan- 
sion for / e C 4+ "(R), where n is the number of error terms, as well as stability 
for the Crank-Nicolson scheme using Gershgorin's theorem to determine the spec- 
trum of the Griinwald matrix. Richardson extrapolation then gave second order 
convergence in space and time of the numerical scheme. 

In this article we explore convergence with error estimates for higher order 
Griinwald-typc approximations of semigroups generated first by a fractional de- 
rivative operator on L (M) and then, using a transference principle, by fractional 
powers of group or semigroup generators on arbitrary Banach spaces. 

It was already shown in [5J Proposition 4.9] that for all 

/ G X a (R) := {/ G Li(R) : 3g G Li(K) with g(k) = {-ik) a f(k),k G M} 
the first order Griinwald scheme 

(1) A?„f(x) = — i— — V T{ - m ~ a ] f(x- (m-p)h) 

y ' m=0 v ; 

converges in to as h 0+ and any shift p G K. Here g(k) = 

e lkx gix) dx denotes the Fourier transform of g and for / G X a (WL), = g iff 
(—ik) a f = g. In Section 3 we develop higher order Griinwald-type approximations 
. In Corollary 13.51 we can then give the consistency error estimate 

(2) A%f-f {a) , <Ch n \\f^\\ Lm 
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for an n-th order scheme. 

Using a new Carlson-type inequality for periodic multipliers developed in Section 
2 (Theorem 12 .4[) we investigate the stability and smoothing of certain approxima- 
tion schemes in Section 4. The main tool is Theorem 14.11 which gives a sufficient 
condition for multipliers associated with difference schemes approximating the frac- 
tional derivative to lead to stable schemes with desirable smoothing. In particular, 
we show in Proposition 14.21 that stability for a numerical scheme using ([TJ to solve 
the Cauchy problem 



(3) 



9 u(t,x) = (-iy +1 -^u(t,x);u(0,x) 



Of 



with 2q — l<a<2q + l,q£N can only be achieved for a unique shift p; i.e. it 
is necessary that p = q for (— 1) 9+1 A% to generate bounded semigroups on Li(K) 
where the bound is uniform in h. Furthermore, in Theorem 14.51 we prove stability 
and smoothing of a second order scheme. 

Developing the theory in L\ allows in Section 5 the transference of the theory 
to fractional powers of the generator —A of a strongly continuous (semi-)group 
G on a Banach space (X, || • ||), noting that f(x — (to — p)h) in (JJ) will read as 
G((m-p)h)f [2]. The abstract Griinwald approximations with the optimal shifts 
generate analytic semigroups, uniformly in h, as shown in Theorem 15.11 This is 
the main property needed in Corollary 15.31 to show that the error between the 
solution S a (t)f — e^ -1 )' A f and a fully discrete solution u n obtained via a 
Runge-Kutta method with stage order s, order r > s + 1 and an JV + 1 order 
Griinwald approximation is bounded by 



\S a (t)f-u n \\ <C 



\.f\\+h 



N+l 



log 



t 



\A N+1 f\\ 



h>0, t = 



Note that this yields error estimates of our numerical approximation schemes ap- 
plied to ([3]) in spaces where the translation semigroup is strongly continuous, such 
as L p (M), 1 < p < oo, BUC(M.), Co(K), etc. Using the abstract setting we can also 
conclude that the consistency error estimate ([2]) holds in those spaces, with the L\ 
norm replaced by the appropriate norm. 

Finally, in Section 6, we give results of some numerical experiments, including a 
third order scheme, highlighting the efficiency of the higher order schemes and the 
dependence of the convergence order on the smoothness of the initial data. 



2. Preliminaries 

A measurable function ip : M. — > C is called an Li-multiplier, if for all / S L± = 
L^R) there exists a v € L\ such that ipf = v where f(k) = e lkx f(x) dx 
denotes the Fourier transform of /. Define an operator, using the uniqueness of 
Fourier transforms, : L\ — > L\ by T^f := v where v is defined as above. It is 
well known that is a closed operator and since it is everywhere defined, by the 
closed graph theorem, it is bounded. Moreover, if ip is an Li-multiplier then ip(k) = 
jj,(k) = e lkx d/i(x) for some bounded Borel measure /i and ||^IIb(Li) = 1 1 A* 1 1 t v 
where TV refers to the total variation of the measure and IHIg^^ is the operator 
norm on L\. Furthermore, if the measure p has a density distribution g then 

(4) \n\\ B[Ll) = 11.911^ = II^IL' 
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where $(x) = w- e~ lkx g(k) dk denotes the inverse Fourier transform of ip. In 
the sequel we will also make use of the following scaling property. If a £ R, h > 0, 
and i[> a (k) = ip(k + a) and iph(k) = ip(hk), then 

( 5 ) II^VIIb(Li) = W^hWeiLt) = W T '4' a lle(Li) • 

The following inequality is of crucial importance in our error analysis. It is a 
special case of a more general Carlson type inequality, see [T31 Theorem 5.10, p. 107]. 
We give an elementary proof here to keep the presentation self contained. The case 
r = 2 is referred to as Carlson-Beurling Inequality, and for a proof see [IJ p. 429] 
and [7]. 

The space W^R), r > 1, denotes the Sobolev space of L r (R)-functions with 
generalised first derivative in L r (R); that is, / G W^R) if / G L r (M.), f is locally 
absolutely continuous, and /' G L r (R). 

Proposition 2.1 (Carlson-type inequality). If g G (R), 1 < r < 2, then there 
exists f £ L\ and a constant C = C(r) > 0, independent of f and g, such that 
f = g ffl.e. and 

(6) ||/|| Ll <C||fl||£>'||i., 
where - + - = 1. 

r s 

Proof. First, since g G L r (R), / := g G L S (R). Also recall the Hausdorff-Young- 
Titchmarsh inequality j6[ p. 211] 

(7) llull, <V2^\\u\\ r ,«eir, l<r<2, - + - = 1. 

r s 

Note that ||/|| Li = ||.g|| Li = ^:||<?|| Ll and 5 '(x) = (-ix)g(x). If g = there 
is nothing to prove. Otherwise, using Holder's inequality, ([7|), and setting v = 

(r-D^n 9 iL r we e et 

ii/ik = ^ iiaik = i( / \m\dx+ 





2tt I 




2 l/r 
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2 l/r 


< 






V2^ 



















-(a$(a:)) 
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-l/r 




(r-1)- 1 / 


') 



2 l / r i i 

_ y 1/r 2 llfflll p llff'IIL • 

Thus / g Li, / exists and g = g = f a.e. □ 
Remark 2.2. Inequality ([6]) can be rewritten in multiplier notation as 

l|r g || B(Ll) <C|| S ||-J| 5 '|lt. 

The reach of Carlson's inequality can be greatly improved by the use of a parti- 
tion of unity: 
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Corollary 2.3. Let <pj be such that J^j 4>j{ x ) = 1 f or almost all x. If gcf>j G W^(M), 

i j_ 

1 < r < 2, for all j and y\ llff^jllf, ll(50j)'ll£ < °°; where i + ^ = 1, i/ien £/iere 
exists f G L\(M.) and a constant C = C(r) > independent of f and g, such that 
f = g a.e. and 



\f\\L 1 <cY J hUi r \\{9^)'\ 



Proof. By design, g = J^j 94>j ■ Let /j G £i(R) be such that fj = g<f)j by Proposition 
12.11 By assumption 



i=i j 



So H'jLifj converges to some / G L 1} < C(r) £V 

/ = E 7 - /j = 9 a - e -: as hmn-^oo /- YJj=i h ^ lim ". 



0. 



u 3\\L r 



"]) \\L, 



and 



□ 



2.1. Periodic multipliers. For periodic multipliers a Carlson type inequality is 
not directly applicable as these are not Fourier transforms of Li-functions. We 
mention that in [5] a suitable smooth cut-off function rj was used where ry = 1 in 
a neighborhood of [— tt, tt] and r\ has compact support to estimate the multiplier 
norm of a periodic multiplier ip by the non-periodic one nip. For the multiplier 
norm of r/tp the above Carlson's type inequality can be then used. However, we 
prove a result similar to Proposition 12.11 for periodic multipliers which makes the 
introduction of a cut-off function superfluous and hence simplifies the technicalities 
in later estimates. 

The space W} per [— tt, tt] below denotes the Sobolev space of 27r-periodic functions 
g on R where both g and its generalized derivative g' belong to L r [— tt, tt]. 

Theorem 2.4. Let g G W} [ — 7r,7r], 1 < r < 2, then g is an L\-multiplier and 
there is C — C(r) > 0, independent of g, such that 



\T g \\ B{Ll ) < \ao\+C\\g\\ 



LA-i 



where - + - = 1 and a = g(x)dx, denotes the th Fourier coefficient of g. 

Proof. Since g G L r [— tt, tt] and is 27r-periodic it can be written as a Fourier series 
9( x ) = J2kL-oo a k elk2 where a k — 57 Jl e~ lkx g(x)dx, denotes the k th Fourier 
coefficient of g. If [i := X)fcL-oo a k5k, where 5k is the Dirac measure at k, then g = jx 
and hence g is an Li-multiplier if and only if || 2"^ || B(z, a ) = IImIItv = Yll^-oo \ a k\ < 
00. First, note that |ao| < \g(x)\dx < 00 and that ika k are the Fourier 

coefficients of g' . Using Bellman's inequality [4] with a = /3 = s, and Hausdorff- 
Young inequality, (see [313 [23]), we have 




\(ika k )\* 



<C\\g\ 



L r \-r, 



\L T \-7, 
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Clearly, the same inequality holds for X)fei-oo l afe l- Thus 

oo 

v— v I 1 

2^ \a k \<\a \ + C ||ff||2 rHr>ir] \\g \\l r[ -* t x] < oo 

k— — oo 

and the proof is complete. □ 

Remark 2.5. The term |ao| cannot be removed from the above estimate in general 
as the specific example g = 1 shows. 



3. Consistency: Higher order Grunwald-type formulae 

Let a £ R+ and let 

X a (R) := {/ G Zi = Li(R) : 3 .g G Li with g(fc) = (-ik) a f(k), k £ R}. 

For / G X Q (R) define /( Q ) = g if g € L x and (-ik) a f(k) = g(k) for fc £ R. 
The function defined uniquely by the uniqueness of the Fourier transform, is 
called the Riemann-Liouville fractional derivative of /. Set ||/|| a := ||/ H^pR) 
for / £ X a (R). Similarly, we define 

X a (R+) := {.f £ Li(R + ) : 3g £ L X (R+) with g(z) = Rez < 0}, 

where f(z) — f^°e zt f(t)dt denotes the Laplace transform of /. For / £ X a (R + ) 
define f^ = g if g £ ii(R+) and (-z) a f(z) = g{z) for Rez < 0. Set j|/|| Q := 

ll/ (a) IL l(R+) for /e^(R + ). 

In order to calculate the fractional derivative of /, a shifted Griinwald formula was 
introduced in |18j given by 

where p is an integer. 

Remark 3.1. When (JS) is applied to a function / G Li(R+), we extend / to Li(R) 
by setting f{x) = for x < 0. Hence, if p < 0, then with this convention, A% f is 
supported on R + and hence can be regarded as an operator on Li(R + ). 

In [21 Proposition 4.9] it is shown that for all / £ X a (R) we have A% p f -> f (a) 
in Li(R) as ft -> 0+. For p = 0, a > and / G X a (R + ), see [HI Theorem 13] for 
the same result. Theorem 13.31 shows that for / £ X a+ p{R) the convergence rate in 
Li(R) is of order /i^,0</3<las/i— s- 0+, and if p < 0, then the same holds for 
/ g X a+P (R+) in Li(R+). 

Furthermore, a detailed error analysis allows for higher order approximations by 
combining Griinwald formulas with different shifts p and accuracy h, cancelling out 
higher order terms. This will be shown in Corollary 13.51 

In the error analysis of Theorem 13.31 below, the function 

(9) oj p>a (z) = 

where a £ R+ and pGl, plays a crucial role. As we take the negative real axis as 
the branch cut for the fractional power, u}p >a is analytic, except where (1 — e~ z )/z 




6 



B. BAEUMER, M. KOVACS, AND H. SANKARANARAYANAN 



is on the negative real axis and at z = 0. As lim z _>.o tu p ^ a (z) = 1, the singularity at 
zero is removable and hence there exists R > 0, a,j £ K such that 

oo 

(10) w P ,a(z) = ^2 a i z ° for all |z| < 

3=0 

In particular, do = 1 and oi = (p — ol/2). Furthermore, this implies that there 
exists C > such that 

f C\z\ for |z| < i?, 

(11) \u Pta {z)-l\< I C forzGiK, 

[ C for Rcz > & p < 0. 

It was shown in [21[ Lemma 2] that 

(12) ^ J =y e- 2 *g Q (t)dt, Rez>0, q Q eii(R+), 

and J °° f/ Q (t) dt = 1. 

Moreover, for a > 1 and fc € R, or for < a < 1 and — ir < k < it, it is easily 
verified that 

(13) 4r(u n .J-ik)) < 
Lemma 3.2. Let 



gp,N, P {k) 



{-ik) N +P 

where w PiQ (z) is given by &j 6j/ (|10[) and 0</3<l,-/V€N. XTien gp,N, P is the 
Fourier transform of some gp,N, P G £i< Furthermore, supp(<7^Ar. p ) C [— p, oo)UR+. 

Proof. Note that for AT, p = 0, /3 < 1, this was shown in [3TJ Lemma 5]. Let N = 0, 
P < 1 and p ^ 0. Then 

&,»,(*) - 1 - + 

The inverse of the first term is given by t i-> g/3 t ofi{t + p); the second term inverts 
to + p) l3 - 1 H(t+p) - t l3 - 1 H(t))/T(/3), where if is the unit step function. 

Clearly, both functions are in L\ and their support is contained in [— p, oo) U R+. 
Assume TV = and /3 = 1, or N > 1. 

Case 1: a > 1 in equation ((9|). Clearly, each term of g is the Fourier transform of a 
tempered distribution satisfying the support condition and hence all that remains to 
show is that gp,N, P is the Fourier transform of an L\ function. Note that gp,N,p tfL L2 
for f3 < 1/2. Hence we use Corollary 12.31 choosing a particular partition of unity 
{<j>j) jeZ . If in ([ID]) R > 1, set R = 1 and define (f> (x) = 1 for \x\ < R/2, 2(R-\x\)/R 
for R/2 < \x\ <R,0 else. Then supp(0 o ) C [-R,R]. Define fafa) = 2{x-R/2)/R 
for R/2 < \x\ < R, <j)x(x) = 1 for R < x < 1 and <j} X {x) = 2 - x for 1 < x < 2 and 
supp(0i) C [R/2, 2}. For j > 1, let 

(x - 2 j ~ 2 )/2i- 2 for 2J- 2 <x< T' 1 , 
(t>j(x) = { (2 J - x)/2J- 1 for 2J- 1 < x < 2*', 
else. 
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Then supp(0j) C [2 J 2 ,2 J ]. For j < 0, define (j>j(x) = <fr-j(—x). Hence there exists 
C such that ||0j-(x)||oo < C2~^+ 2 for all j G Z. By Corollary El we are done if 

we can show that J2- Wd4>j ll£ 2 II (ff^j)'ll£ 2 < 00 ■ 

As .g is analytic in (-2i?, 2i?), <?</> G W%(M). For j 7^ 0, 



90,N, P ( k m( k ) = (-ik)N+p ™ ~ j^ikY^^ 

=Ti(fc)-T 2 (fc). 



By ([TT]) and (|13[) . and cj^ q are bounded. Recall that either N = and 
/? = 1, or N > 1. In either case, the exponent in the denominator of T\ is at 
least one. Hence there exists C independent of j such that |Xi(fc)| < C2 2 ~> : >>, and 
\T[(k)\ < C2 2 -I J 'I. The length of the support of T x is less than 2^ . Hence 



|Ti||| < c 1 / 2 (2 4 - 2 l J l2l J l) 1 / 4 < C2-I J 



I /-1 



Same for \\T{\\1 2 . To estimate recall that ||<^-||oo < 2~W+ 2 . Thus there exists 

C independent of j such that \T£(k)\ < C2^\- 2 ^-^ + C2^- 2 ^-^2-^+ 2 = 
C2 3+p 2 -(p+i)\j\ and hence 



II^IIL < C(2- 2 ^ +1 ^2^) 1 / 4 = C2-\^ 2f}+1 ^\ 
The norm of Ti is bounded by 

||T 2 ||1 < C(2- 2 l J ^2l J l) 1 / 4 < C^ 1 - 2 ^ 4 . 



This implies that \\g^N,p4> 3 \\L 2 \\(gp,N, P ^y\\L 2 < C2\^~ 2 ^l 4 2~\^ 4 = C2-^/ 2 , 
and hence J2j \\gp,N,p4>j\\L 2 \\{gp,N,p4>j)'\\L 2 < 00 and therefore gp,N, P G Li(M). 

Case 2: a < 1 in equation ©. For a < 1, w^, a is not bounded; if a < 1/2 
it is not even locally in L2, making the above method unfeasible. Instead we use 
induction on N. We already established that the assertion holds for N = 0, /? < 1, 
so assume it holds for some TV; i.e. assume <jp,N,p is the Fourier transform of an L\ 
function satisfying the support condition. Then using the fact that the convolution 
of L\ functions is an L\ function and the fact that we established the assertion for 

a = 1, we obtain that k 1— > gp,N,p{k) — — is also the Fourier of an L\ function 
satisfying the support condition as the support of the inverse of the second factor 
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is in R+. Furthermore, 



9p,N, P (k) 



—ik 



- 1 



—ik 



-ikp 



—ik 



-ik 



(-^fc) JV + 1 +' 3 
(^) a e-^-Ef=o^H^ 



( l-e ik 
\ -ik 



{-lk) N + l +f> 



a 

' 1-e 



-iA- 



(_^)A+l+/3 



--h-h-h, 



where the are Taylor coefficients of w p , Q +i, the c m are the Taylor coefficients of 
— | — and d, are such that equality holds. As by case 1, I\ and I2 are Fourier 
transforms of L\ functions, so is I3. Hence ^3 at k = has to be bounded and 
therefore dj = a,j and hence gp t N+i,p = k- 

The same argument about the convolution of Li-functions applies to the remain- 
ing case of N = 0, (3 = 1, using (fl"2|) and the fact that 



U)p, a (-ik)- 



-ik 



1 



■>p,a+i(-ik) - 1 w p<a (-ik)-l 



-ik 



-ik 



-ik 



□ 



Theorem 3.3. Let < j3 < 1. TTien i/iere exists C > suc/i i/ia£ / e ^ a+ ^( 
implies that 



Ah, P f f 



(a) 



ill 



<C7/||/IU, 



ash ^ 0+. 7/p < and / e A" Q+( g(R + ), fftera 



^, P / - / (a) 



<C^||/| 



/n case of /3 = 0, / £ X Q . 



< P / - / (a) 



Li(K) 



0. 



Proof. The case j3 = was already shown in [2 Proposition 4.9] and is only included 
in the statement of the theorem for completeness. 

To prove the first statement, first note that (1 + z) a = X)m=o (m) z?n ^ or z e 
C, \z\ < 1 and a € M.+ . Moreover, the binomial coefficients are related to the 
Gamma function by the following equation = r(- n )r(m+ij ■ Taking Fourier 
transforms in ([8]) we get 

{AfJ){k) =h- a J2 {-l) m [ a )e %k{m - p)h f{k) = h- a e- lkhp {l - e lkh ) a f{k) 

m=0 \ m / 

={-ik) a u v , a {-ikh)f{k) = + 
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where uj PtCC (z) is given by ((9|) and (h(k) = (—ik) a (ujp ta (~ikh) — l)f(k). We rewrite 
the error term Ch(k) as a product of Fourier transforms of Li-functions: 

(14) C h {k) = (-ik) a (u: p , a (-ikh) - l)f(k) = h^gp, 0lP (kh)f^+P)(k), 

where, by Lemma 13.21 gpfl. p {k) = (^ ulp '^ i ^) — -j is the Fourier transform of an 

Li function gp,o,p and /( Q +' 9 ) g L\ by assumption. Thus is indeed the Fourier 
transform of an Li-function ^ and, by Q, ([5]) and Lemma [ 



|0i|| Ll (R) < h P \\T-g 0iOiV {h.)\\B( Ll { 
<Ch^\\f\\ a+fj . 



>+/3) 



hP \\9Pfi,v\\ Ll ( 



(«+/3) 



Finally, the second statement follows the same way, taking Remark 13 . 1 1 into account 
and noting that if p < 0, then supp(£/i) C R+ and using Lemma T3. 2 1 in (|14[) . □ 

Remark 3.4. Let a G R + , 2g — 1 < a < 2q + 1, ?eM, and 
(15) 1> a ,hAz) = (-l) q+1 h- a e- h P z (l - e hz ) a . 

Then the proof of Theorem 13.31 shows that the Griinwald formula can be ex- 
pressed in the multiplier notation of Scction[5]as A% p = T^_-^ q +i^p where 



(16) 



^ h (k) = iP a ,h, P m = (-l) q+1 h- a e- lkh P(l - e lkh Y 



Analyzing the error term in (|14[) further allows for higher order approxi- 
mations by combining Griinwald formulae with different shifts p and accuracy h, 
cancelling out lower order error terms. Consider the Taylor expansion of w PjQ , in 
(flOl and let N > and for < j < N let bj,Pj G R and Cj > be such that there 



exist dj with 



N 



3=0 



j=N+l 



for \z\ < 2R; i.e. consider a linear combinations of w PjiQ ,'s cancelling out the lower 
order terms. Define 



(17) 



3=0 



to be anW + 1 order Griinwald approximation. This is justified according to the 
following Corollary. 

Corollary 3.5. Let < /? < 1, TV G N and Jl£ 6e an iV + 1 order Griinwald 
approximation. Then there exists C > such that / G X Q+ 7v+/i(R) implies that 



(a) 



/HcH 



as h ->• 0+. < /or all < j < N and f G X Q+ 7v+/3(R+) i/ien 



(n) 



< C/i^ll/IUw 
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Proof. Note that ^ = anc ^ f° r > 0, we have that ^2j =0 bjaj. n dj = 0. 

Following the same argument that led to (fT4)) . we obtain that the error term can 
be expressed as 

N 



Uk) =(Alf - /(«)) (k) =X) 6 J K.aHAcj-fc) - 1) (-ik) a f(k) 

AT / TV \ 

= ^6j u Pj>a (-ikcjh) - ^ a,j,n(-ikcjh) n J (-ik) a f{k) 

3=0 V n=0 / 

j=0 v J y 



where are the Taylor coefficients of ui Pj ^ a . By Lemma 13.21 



is the finite sum of Fourier transforms of L\ functions and hence is the Fourier 
transform of an Li function with 

llalU l( M)<C^||/|U + Ar +/3 . 

The second statement follows along the same lines. □ 

4. Stability and Smoothing: Semigroups generated by periodic 
multipliers approximating the fractional derivative 

The next result is the main technical tool of the paper. It gives a sufficient condi- 
tion for multipliers associated with difference schemes approximating the fractional 
derivative to lead to stable schemes with desirable smoothing. In error estimates 
later, the smoothing of the schemes will be used in an essential way to reduce 
the regularity requirements on the initial data to obtain optimal convergence rates 
when considering space-time discretizations of Cauchy problems with fractional 
derivatives or, more generally, fractional powers of operators. 

Theorem 4.1. Let a £ R+ and tj) be an absolutely continuous 2tt -periodic function 
that satisfies the following: 

(i) < C\k\ a for some C> 0, 

(ii) \i>'(k)\ < C'\k\ a ~ X for someC > 0, 

(iii) Re('0(fc)) < — c \k\ a for some c > 0. 

Then t/j £ W^ per [— it, tt], where r = 2 if a > ^ and r < — if a < |. Moreover, 

(a) ||T e *|| B(il) <Kfort>0, 

(b) < f fort >0, 

where K and M depend on c, C and C above. 

Proof. Set r = 2 if a > | and r < j^- if a < | , where 

(18) 1 < r < 2 and - + - = 1. 

r s 
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Note that r(a — 1) > —1, then from the assumptions we have ip £ W^ per [— 7r, it], 
and for t > 0, we also have G H^ 1 [ — 7T, 7r]. By Theorem 12. 4i 



(19) 



llr^ll^^laol + Clle^H^IKe^yi^ 

Firstly |ao| - \^J^e^ k Uk\ < J-j^e tR <^dk < & JZ„ e^M" dk < 1, 
where we have used Assumption (iii). Now, using Assumption (iii) again together 
with the substitution r = t~\k\, 



lL r [-7r,7rl 



dk 

-rct\k\ a Ak 



< C 



s rtRe(V(fe)) dfc 



(20) <C 

Making use of Assumptions (ii) and (iii) , we have 



A. ( e t^)\ 


r 




' < C (t r |fcf (a-1) ) 


dk \ ) 




dk 





Since r (a — 1) > —1, an application of (fT5|) and the substitution r = £<*|fe|, yields 



(21) 



<c t 



L r \ — 7T,7r| 



r r(a-l) e -rcr° dr \ < ^ 



The proof of (a) is complete in view of (JT^, (J5DJ) and ([2"Tj). 

We have that ^e'^ G W^ per [—TT, n] and once again by Theorem l2.4 



(22) 



Note that |o | < ^ |V>(*0I e tRe W fc » dfc < C* |fc| Q e ~ c *l fc r dfc < C/t, where 
we have used Assumptions (i) and (iii). The use of the substitution r = t«\k\ and 
the Assumptions (i) and (iii), yield 



(23) 



\ipe 



ti/i\ 



\LA-n 



<C[ |fc|™ e - rct|fc| ° dfc ) " < Ct 



We also have by virtue of (|2ip and the three assumptions, 



dk 



(V^Xfc) 



dj>{k) 
dk 



<2 r - i (t r \ip(k)\ r +i)\(my 



< C{\k\ r{a - 1) +t r \k\ r(2a ~ 1) )e- rcm " . 

Thus, using (fT5)l and the substitution r = ti |fc|, and noting that r(2a — 1) > 
r(a - 1) > — 1, 



(| r |K«-i) + | r r( 2Q - 1 ))e- rcTQ dT 



L r [— 7T,7r] 

(24) <C*t-^+^, 

and the proof of (b) is complete in view of (|T8]), (|22|) . ([23]l and ([24 



□ 
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4.1. The shifted Griinwald formula of order 1. First, we make the following 
two observations about the multiplier associated with the shifted Griinwald formula 
given by (fT6"|). Note that for p, k e K and h,a £ R + , such that 2q — 1 < a < 
25+I, q€N, 

(25) Vh(*0 = (-l) 9+1 ^ Q e-' fc ' ip (l - e 4fc,l ) a = (-l) q+1 h- a (-ikh) a oj p>a (-ikh), 
where w Pi0 ,(2) is given by ©. Note also that 

(26) ip p h (k) = h- a i}[{hk). 

We now show that the range of the symbol associated with the shifted Griinwald 
formula is completely contained in a half-plane if and only if the shift is optimal. 

Proposition 4.2. Leti> p h (k) = (-l)q+ 1 h- a e- iphk (l-e ikh ) a , where the shift p£ N, 
k e R, h > 0, a € R+ and q E N such that 2q— 1 < a < 2g+ 1. Then 

(a) ^ satisfies Assumptions (i) and (ii) of Theorem \4-l\ with C, C indepen- 
dent of h. 

(b) Re(^) does not change sign if and only if |p — |H < ^ «/ and only if ip^ 
satisfies the Assumption (Hi) of Theorem \4- 1\ with c independent of h. 

Proof. Proof of (a) : 

Recall that by dTTJ and ([T3|l. \u Pt0l (-ikh)\ < C and \aj' pa (-ikh)\ < C, for 
k E [— vr, 7r], respectively. Thus using ([23)1 . for fc 6 [— 7r,7r], 

(27) |^(fc)|<C|fcr,and 
(28) 

-aih{-ikh) a - x u Pta {-ikh) + {-ikh) a u' pa {-ikh)\ < C^ 01 ' 1 , 



dr h (k) 



dk 

for some C > 0. 

Proof of (b): It is enough to consider the case h = 1 in view of (|26|) . so let 
V' := V->i- Since ^ is 27r-pcriodic and tp(k) = ip{—k), it is sufficient to consider 
k e [0,tt]. Now 

^(fc) = (-l)«+ 1 e-** (e*(e^ - e T))° = (-l)«+ 1 e < « ~ p)fc f-2isin . 
Using the fact that for x > 0, (— ix) a — x a e _M f , we obtain 

el ( f _p) fc _l-L 0<A , <7ri 



^(fc) = (-l) 9+1 2 Q sin« 



and therefore 



Re (m) =(-l) 9+1 2 Q sin" - cos ((£ - p)fc 

(29) 



OTT 



=( _ 1)9+ i-r 2 a sin . / k \ _ _ ?r) 



Clearly, as < fc < 7T, (f2"9")l changes sign if and only if | ^ — p > |. Note that by 
assumption |^— p| 7^ 5. Furthermore, Assumption (hi) of Theorem 14. II implies that 
there is no sign change. Hence all that remains to show is that \p — ^ | < | implies 
Assumption (iii) of Theorem 14.11 The fact that c is independent of h follows from 
|J2BJ|. 
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Note that \p — ^ | < \ implies that p = q and hence, using that for < x < it, 
s'm(x/2) > x/ir and cos((% — p)(x — it)) > cos(— {% ~ p)ir) , 



Re (m) = - 2 Q sin" ( - I cos [{- - p)(k - tt] 



a 



(30) / k 

<-2 a cos ((| - p)n) = -k a 2 a cos ((| - P )tt) /tt". 



□ 



Next, we show that with the optimal shift; i.e., p = q, the operators T^ h generate 
strongly continuous semigroups on Li(R), and in the case whenp = 0; i.e., < a < 
1, on that are bounded uniformly in ft. Since the range of iph is always 

contained in the spectrum of T^ h , the semigroups generated by T^ h will not be 
uniformly bounded if the shift is not optimal as shown by Proposition 14.21 In fact 
we show more: if the shift is optimal then the semigroups are uniformly analytic 
in ft; i.e., there is M > such that the uniform estimate \\T^, he t^ h \\b(Li) < Mt~ x 
holds for t,h > 0. This fact will have significance when proving error estimates for 
numerical schemes for fractional differential equations. 

Theorem 4.3. Let a eW. + , 2p — 1 < a < 2p + 1, p G N, and 

ip h (k) = (-l) p+1 h- a e- ipkh (1 - e lkh ) a . 

Then the following hold. 

(a) {T e t^ h }t>o are strongly continuous semigroups on Li(R) that are bounded 
uniformly in h > and t > 0. In particular, if 1 < a < 2, then {T e ti, h }t>o 
is a positive contraction semigroup on Li(R) and for < a < 1, on Li(R + ). 

(b) The semigroups {T e t^ h } t >o are uniformly analytic in h > 0; i.e. there is 
M > such that the uniform estimate \\T^ he t^ h \\b{l±) ^ Mt holds for 
t,h>0. 

Proof. Proof of (a): To begin, note that iph(k) = i>h(k) with p = q. We only have 
to show that ||T e t^ h Herz-i) — ^ or au ^ — ^ an< ^ h> Q, for some K > 1 and strong 
continuity follows by [TJ Proposition 8.1.3]. Furthermore, it is enough to consider 
ft = 1 in view of ([5]) and ([25]), so let tp := ipi- First, lct0<a<lorl<a<2 
and hence p = or 1, respectively. We have, taking Remark 13.41 into account, that 



(T^f){x) = (_l)P+i(f;(_l)™["V( s _ (m-p))) 

m=0 ^ ' 

l)m+(-iy +1 E (-ir(° 



f(x - (m-p)) 



a 



= (-{)lf)(x) + (Tjf)(x). 

Since (— l) p+1 (— l) m ( r ") > 0, for m ^ p, it follows that Tj, is a positive operator on 
Li(M) (or, Li(R+) for < a < 1 by recalling Remark l3~Tj) and so is e tT * — T etn j,. 
Therefore, noting the fact that ^ =0 (-l) ra Q = 0, 

rp tT^ t(~( a )l+T 7 ) -( a )t tT 7 ^ n 
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and 

II^Hl, < e -G)*e*ll T *"^ = e -C)* e (- 1 ) P+1 ^=o, m ^(-i) m (™)* = i. 

Let now a > 2, then by Proposition 221 i>h satisfies the hypothesis of Theorem 14. II 
and the proof of (a) is complete. 

Proof of (b): Let a > 0, then the statement follows from Theorem 14. II in view of 
Proposition 14.21 

□ 

4.2. Examples of second order stable Griinwald-type formulae. Let a £ 

R_l_, 2q — 1 < a < 2<jf + 1, and q € N. Consider the mixture of symbols of Griinwald 
formulae yielding a second order approximation of the form 

4> h {k) + 

where ij) p h (k) = (-l)i +1 h~ a e- ipkh (l - e lkh ) a , the symbol of the p-shifted Griinwald 
formula and thus = (—l) q+1 T ( f >h . There are of course many combinations of 
a,a,pi, and pi that yield a second order approximation; however only some are 
stable. It is straight-forward to show that if < a < 1, then a = 2,pi = P2 =0 
give a second order approximation and similarly for l<a<2,a = 2— — , p± = 
lj P2 — \- So we only show stability. 

Proposition 4.4. Let <f>h(k) be as above, where a — 2, p\ = p2 = 0, if < a < 1 

and a = 2 — ^, pi = 1, P2 = if 1 < a < 2. Then <frh satisfies the assumptions of 
Theorem \4-l \ with constants c, C and C independent of h. 

Proof. Assumptions (i) and (ii) of Theorem 14. II are clearly satisfied in view of (|2"T)) 
and (j28j) . 

Proof of (iii): As 

Mk) = h- a (/ )1 {hk), 
it is sufficient to consider the case h = 1, so let <j> := <f>\. That is, 

<j>{k) = 0i (k) = (ae- tpik {l - e tk ) a + (1 - a)2~ a e~ t2p2k {\ - e t2k ) a ) 

and note that <f> is 27r-periodic. By symmetry, we only need to consider the real 
part for < k < ir, and so using (|29| and the double angle formula we have 

Rc(0(fc)) = (-l) 9+1 2 Q sin a ^acos^+(l-a)cos" ^ cosS 

where A = (f - Pl )k - Sf and B = (a - 2p 2 )k - *f. 
We consider for < k < ir, the function 

F(k) = {-l) q+1 {acosA+ (1 -a)cosB) 

and show that F(k) < F(ir) < 0. Then, since a > and p\ is the optimal 
shift, if Re ((1 - a)V>f (fe)) < 0, by (|30]), Re(0(ft)) < -ak a 2 a cos ((§ - p)n) /tt q . If 
Re((l — a)ip P ^{k)) > 0, we will have the estimate 

Re(0(fc)) < F(7r)2 Q sin Q (*) < k a F '{ir)2 a / 'ir a . 
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An easy check shows that F(tt) < 0. It remains to show that F'(k) > 0. Now, 
as in both cases 2pi = pi and (1 — o)(q — 2p2) = —a — p±), 

F\k) =(-l) 9+1 (-a (| - Pl ) sin(A) - (1 - a) (a - 2p 2 ) sm{B)) 
=C(a) (sin(A) - sin(B)) = 2C{a) cos 

where C{a) = -a(f - Pl ), A±b = - pi)k - eg, and A^R = Checking 
the range of the arguments, the cosine factor is negative if a > 1 and positive if 
a < 1, the sine factor is always negative, and hence F'(k) > 0. □ 

Theorem 4.5. Let <f>h{k) and a be as in Proposition \4-4\ If < a < I and a = 2, 

or if 1 < a < 2 and a — 2 — — , then {T e t* h } t >o are semigroups on Li(R_|_), or 
ii(M), respectively, that are uniformly bounded in h and t, strongly continuous, 
and uniformly analytic in h. 

Proof. The statement follows from Theorem 14. II in view of Proposition POl □ 

5. Application to fractional powers of operators 

Let X be a Banach space and —A be the generator of a strongly continuous 
group of bounded linear operators {G(£)}t 6 R on X with ||G(i)||g(x) < M for all 
t G R for some M > 1. If fj, is a bounded Borel measure on M and if we set 
-0(z) := = e zs d/i(s), (z = ik), we may define the bounded linear operator 

(31) ip{-A)x:= I G(s)xdn(s), x e X. 

It is well known that the map ip ip(—A) is an algebra homomorphism and is 
called the Hillc-Phillips functional calculus, see for example, [11]. That is, if ip = fi 
and (f) = z), for some bounded Borel measures /z and v, then (</> + ip){— A) = 
4>{-A) + ijj{-A), {(j> ■ tp)(-A) = c/)(—A)tp(—A) and {c<p)(-A) = c<f>(-A), c 6 C. A 
simple transference principle shows, see e.g. [21 Theorem 3.1], that, 

(32) II^M)||b(jc) < M||T fe ^ (ifc) || B(il(R)) . 

Note that if supp/z c K + , then we may take — A to be the generator of a strongly 
continuous semigroup and the properties of the Hillc-Phillips functional calculus 
([31]) and the transference principle (|32[) still holds. 

Let 2p — 1 < a < 2p + 1, a £ R+, peN and {nt}t>o be the family of Borel 
measures on R such that fit{z) = e t (~ 1 * >F+ (~ z ) , (z = ik). Then the operator family 
given by 

S a (t)x:= J G(s)xdn t (s), x&X,t>0, 

is a uniformly bounded (analytic) semigroup of bounded linear operators on X, 
see [U Theorems 4.1 and 4.6] for the group case. In case < a < 1, we have 
slippy C M + and hence — A is allowed to be a semigroup generator and G to be a 
strongly continuous semigroup and the analyticity of S a holds, see [3] and [22] . The 
fractional power A a of A is then defined to be the generator of S a multiplied by 
(— l) p+1 . We note that the fractional power of A may be defined via an unbounded 
functional calculus for group generators (or, semigroup generators), formally given 
by f a (—A), where f a (z) = (—z) a . This coincides with the definition given here 
for groups and, in case < a < 1, for semigroups, see [2] and [3] for more details. 



A + B\ . A-B\ 
[ — ) Sm I —J 
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Thus, for the additional case of —A being a semigroup generator and a > 1, we 
just set A a = f a (—A) as in [3]. 

The following theorem shows the rate of convergence for the Grunwald formula 
approximating fractional powers of operators in this general setting. For the sake 
of notational simplicity we only give the first order version; the higher order version 
follows exactly along the same lines and is discussed in Corollary 15.31 



Theorem 5.1. Let X be a Banach space and p e N. Assume that —A is the 
generator of a strongly continuous group, in case p — 0, semigroup, of uniformly 
bounded linear operators {G(t)} t6 u on X . Define 

Kh x = h ~° E * G « m ~ P) h > = (-^) q+1 ^j l . P (-A)x, xeX, 

m=0 



where ip a h,p{z) * s given by (|15p and a G K + such that 2q — 1 < a < 2<j> + 1, q G N. 
Then, as h > 0, we have 

(33) W a ,h x ~ A ° X W ^ Ch\\A a+1 x\\, x g V{A a+1 ). 

Furthermore, if p = q, then (— l) p+1 $^, h generate {S^ ^(i)}t>o? strongly continuous 
semigroups of linear operators on X that are uniformly bounded in h > and t > 
and uniformly analytic in h > 0; i.e., there is M > such that \\S^ /l (i)||e(A) < M 
and WKAhWUm < Mi" 1 for all t,h>0. 

Proof. If g is defined by g(z) = Up - a ^_^ — i ; Re z < 0, then by Lemma [3.21 with 
N = and (3 = 1 and © we have that \\T kl ^ h ^ ikh) \\ B ^ Ll{R ^ < Ch. In case p = 0, 
we have that supp(.g) C K+ and hence ||T fch 4. h5 ( ife h) ||b(Li(r+)) < Ch. Therefore, 
by the transference estimate (|3"2"j) . /ij4)||g(^) < Cft for some C > 0. Thus, 

if x G X>(A Q+1 ), then ||/i5(-/iA)A q+1 .t|| < C7i||A a+1 a;||. Using the unbounded 
functional calculus developed in [2] (in casep = see [3]), we have, for x G T>(A a+1 ), 



hg(-hA)A a+1 x = hg{hz)( 



\a+l\ 



,'-1 V" I " \ p( m_ P)' l ( z ) 

m=0 



/// 



( — i)"M \ e (m ~ p ' (z) — (—z) a 



z=-A-i 



*l h x-A«x 



and the proof of is complete. 

The strong continuity of h follow from Theorem 14.31 and [H Theorem 4.1], 
where the latter theorem establishes the transference of strong continuity, from 
Li(R) to a general Banach space X (see, [3j Theorem 5.1] for the same result 
in the unilateral case). Finally, the operator norm estimates follow from the -Li- 
norm estimates in Theorem 14.31 and the transference estimate (|32[) , noting that 
by the functional calculus of [2] and [3] it follows that (ip a .h, P e t ^' a - h - p ){— A) = 
(-l) p+1 $ p ath Sl h (t) for t > where ip atKp is given by and p = q. □ 

The stability and consistency estimates of Theorem 15.11 allow us to obtain un- 
conditionally convergent numerical schemes for the associated Cauchy problem in 
the abstract setting together with error estimates. To demonstrate this, we use the 
optimally shifted first order Grunwald scheme as "spatial" approximation together 
with a first order scheme for time stepping, the Backward (Implicit) Euler scheme, 
to match the spatial order. Let 2p — 1 < a < 2p + 1, a G K + , p G N, let X be a 
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Banach space and — A be the generator of a uniformly bounded strongly continuous 
group (semigroup if p = 0) of operators on X and set A a := (— l) p+1 A a . Consider 
the abstract Cauchy problem 

u(t) = A a u(t); u(0) = x, 

with solution operator {S a (t)} t >o, where, as we already mentioned, S a is a uni- 
formly bounded analytic semigroup as shown in [5J Theorem 4.6] and in |22j for 
< a < 1; that is when —A is a semigroup generator. For its numerical approxi- 
mation set 

= (-1) $a h u n+r, ua = x, n = 0, 1, 2, 

T 

that is, with A a , h := {-l) p+1 ^ h , 

u n = (I - TA ath )~ n x, n = 1, 2, ... 
We have the following smooth data error estimate. 

Theorem 5.2. Let 2p - 1 < a < 2p + 1, a G R+, peN, n e N, and < e < 1. 
Let X be a Banach space and —A be the generator of a uniformly bounded strongly 
continuous group (semigroup if p = 0) of operators on X and set t = nr. If 
x S V{A 1+£ ), then 

(34) \\S a {t)x-u n \\ < C^llxlj +h— \\A 1+e x\\), ra = l,2,...;f >0, 

£ 

and, if x G V(A), then 

t 



(35) \\S a (t)x - u n \\ < C(n _1 ||x|| + (1 + a)h 



log , 

h a 



\\Ax\\), n = 1,2, ...;t> 0. 



Proof. To show (|34|) . we split the error as 

S a (t)a; - u„ = S a (t)x - S a ,h(t)x + S p M {t)x - u n := ei + e 2 . 

It was shown in Theorem 15.11 that S p h are bounded analytic semigroups on X, 
uniformly in h. Therefore, ||e2|| < Cn _1 ||x||, n G N, as shown in [8], with C 
independent of h and t. To bound e\ we use the fact that all operators appearing 
commute being functions of A, to write 

(36) ei = S (t)z - S p h (t) = f (A a - A a . h )S a (r)S p h (t - r)xdr 

Jo 

Note that the analyticity of the semigroup S a implies that there is a constant M 
such that for < e < 1, the estimate \\Al~ E S a {t)\\ < Alt 6 ' 1 holds for all t > 0. 
Then, by Theorem 15. 1[ 

\\ei\\<Ch / '\\A a+1 S a (r)S^ h (t-r)x\\dr 
Jo 

= Ch [ \\A^*S a (r)£% h (t-r)A 1 +*x\\dr<Ch—\\A 1+e xl 
Jo £ 

which completes the proof of (|34p . 
To show (f3"5j) , write e\ in (j3"6"j) as 

ei= / (^a-Ax,fc)5„(r)SS h (t-r)a;(ir+ / (4, - A»,h)S a (r)^ h (* - r)xdr 
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It is already known that A a< hX — > A a x as h — > 0+ for all .t <G 2?(A q ), see [2j 
Proposition 4.9] and [21], and hence we have stability ||A Q ./ l a; — A a x\\ < C\\A a x\\ 
for all x £ T>(A a ). Therefore, 



kill <c 

<C{ha 



\A l ~ 



log 



: S a (r)Ax\\dr 
t 



Ch 



\A a S a (r)Ax\\dr 



h° 



\Ax\\ 



□ 



Note that the condition x £ T>(A 1+E ) in might be hard to check for e ^ 1, 
depending on A and the Banach space X . However, one can always use e = 1 as 
T>(A 2 ) is usually quite explicit. We also obtain convergence and error estimates 
of stable higher order schemes (such as the second order Griinwald formulae intro- 
duced in Section l4~2l . 



Corollary 5.3. Let a £ 



with 2q — l<a<2q + l, q£N and let 

N N 
3=0 j=Q 

be an N + 1-order Griinwald approximation, where ip a .h,p( z ) is given by (|15[) and 
bj, Cj,pj are as defined in (|17p . Assume the multiplier X^^Lo ^j^c 3 ?i(^)> w ^ lere "tphik) 
is given by (fTS)) , satisfies (i)-(iii) of Theorem VS. 3\ with constants independent of h. 
If one solves the Cauchy problem 

u{t) = $ a , h u(t); u(0) = x, 

with a strongly A-stable Runge-Kutta method of with stage order s and order r > 
s + 1, then, denoting the discrete solution by u n at time level t = nr, 

t 



\S a (t)x 



< C 



log 



h° 



\A 



N+l r 



|| , h > 0, t 



for all x £ V(A N+1 ). 

Proof. It is straight forward to see that Theorem 15. II holds for ^a^; i.e., 
\\* a , h x ~ A a x\\ < C/i JV+1 ||A 0,+Ar+1 a:|| J x £ V(A a+N+1 ), 



< M and ||* a ,he t * Q ' l || B ( X) < M/t. Following the proof of Theorem 
we obtain the "spatial" error estimate 



(37) 



\\S a (t)x - e tv "' h x\\ < Ch 



N+l 



log 



t 



\\A N+1 x\\ 



Since the analyticity of the semigroups e'* Q ' h is uniform in h (the constant M does 
not depend on h), the statement follows from JTSJ Theorem 3.2] (see also [33]). □ 



Remark 5.4. The error estimates (|34[) and (|35[) arc almost optimal in terms of the 
regularity of the data. We conjecture that one could remove the slight growth in 
t from (|34p or the logarithmic factor in (|35j) by considering the Li-case again and 
using the theory of Fourier multipliers on Besov spaces. Then use the transference 
principle to derive the abstract result. We do not pursue this issue here any further. 
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Remark 5.5. The convergence rate given in Corollary 15.31 can be extended using 
the stability estimate 

\\S a (t)x - e^^xW <C\\x\\ 

and (|37]l to certain real interpolation spaces as in [HI Corollary 4.4]. We note 
that while the spaces V{A S ) endowed with the graph norm are, in general, not 
interpolation spaces, they are embedded within appropriate interpolation spaces 
(see, for example, [5J Corollary 6.6.3]) and therefore we obtain 



\S a (t)x - U n \\ < C 



log 



h° 




, h > 0, t = TIT, 



for all x <G r D(A s ) 1 s G [0, A + 1]. Also note that we indeed have convergence of u n 
to S a (t)x for all x £E X as t — V and h — > by Lax's Equivalence Theorem as a 
consequence of stability and consistency. The order of convergence, however, might 
be very low depending on x. 



6. Numerical Experiments 

In this section we give the results of two numerical experiments. The first is to 
explore the effect of the regularity of the initial distribution on the rate of conver- 
gence as needed for Corollary 15.31 the other is to see how well a second and third 
order scheme fare in the numerical experiment done by Tadjeran et. al |20j . 




10 1 10 2 10 3 



1/h 



Figure 1. Li-error for different initial conditions /j and a first 
(•) and second order (*) scheme. Note the less than first order 
convergence for a "bad" initial condition; i.e. one that is not in 
the domain of A. Also note the less than second order convergence 
for a second order scheme but first order convergence for the first 
order scheme for an initial condition that is in the domain of A but 
not in the domain of A 2 . 

Example 6.1. We consider X ~ £i[0, 1] and A = (d/dx) with 

D(A) = {/:/'eL 1 ,/(0) = 0} 
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and a = 0.8. Wc approximate the solution to the Cauchy problems 
u'(t) = -A a u(t);u(0) = f u i =1,2,3 

at t = 1 with 

h{x)=x-°-\ f 2 (x)=x - 7 , f 3 (x)=x 1 - 7 , 

with first and second order Griinwald schemes (as in Proposition I4.4[) as well as 
via a convolution of fi with an a-stable density, which gives the exact solution but 
both the convolution and the density are computed numerically on a very fine grid. 
Note that f x D(A), f 2 G D(A) but f 2 & D(A 2 ) and f 3 G D(A 2 ). However, 
/i G D(A Pl ) for /3i < 0.7 and f 2 G D(A f32 ) for /3 2 < 1.7. By Remark El] we expect 
about 0.7-ordcr convergence for both schemes in case of uo = fi, and first order 
convergence for the first order scheme for the other initial conditions. We expect 
about 1.7-order convergence for the second order scheme in case of u = f 2 and 
second order convergence in case of uq = f$. For the temporal discretization we use 
MATLAB's ode45, a fourth-order Runge-Kutta method with a forced high degree 
of accuracy in order to investigate the pure spatial discretization error. We see in 
Figure [1] that we obtain the expected convergence in all cases. 

Example 6.2. Even though our theoretical framework is not directly applicable, 
because the fractional differential operator appearing in (|38|) is defined on a finite 
domain with boundary conditions and has a multiplicative perturbation and hence 
it is not a fractional power of an auxiliary operator, we apply the second and third 
order approximations to the problem investigated by Tadjeran et al. [20] . namely 
approximating the solution to 

du(x,t) r(2.2) a , 8 fl L V»,t) _ t 3 , 3 

( 38 ) dt = -^]~ x dx i.8 (l+a)e x-u{x,0)=x 

on the interval [0, 1] with boundary conditions u(Q,t) = 0, u(l,t) = e _t . The exact 
solution is given by e~ t x 3 , which can be verified directly. 

A second order approximation of the fractional derivative is given by Proposition 
14.41 In order to obtain a third order approximation we consider 

<t> h {k) =aip 1 h + inPl h + cip° h 

with the coefficients a, b and c such that <fih is a third order approximation; i.e. 

7-8a + 3a 2 , -7 + 3a 
3(a — 1) 3(a — 1) 

A quick plot of 4>h(k) for k G R strengthens the conjecture that, for a = 1.8, the 
spectrum is in a sector in the left half plane and hence we expect stability and 
smoothing. 

We use again a fourth order Runge-Kutta method to solve the systems to t = 1 . 
Table 1 suggests that we indeed have second and third order convergence with 
respect to the spatial discretization parameter Ax. 

References 

1. W. Arendt, C. J. K. Batty, M. Hieber, and F. Neubrander, Vector-valued Laplace transforms 
and Cauchy problems, Birkhauscr Verlag, 2000. 

2. B. Baeumer, M. Haasc amd M. Kovacs, Unbounded functional calculus for bounded groups 
with applications, Journal of Evolution Equations 9 (2009), no. 1, 171-195. 



HIGHER ORDER GRUNWALD APPROXIMATIONS 21 
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